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SUMMARY 


In  analysing  the  orbit  of  Ariel  1 to  determine  upper-atmosphere  winds,  it 
was  observed  that  the  orbital  inclination  underwent  a noticeable  perturbation 
in  November  1969  at  the  29:2  resonance  with  the  Earth's  gravitational  field, 
when  the  satellite  track  over  the  Earth  repeats  every  2 days  after  29  revolu- 
tions. The  variations  in  the  inclination  and  eccentricity  of  the  orbit  between 
July  1969  and  February  1970  have  now  been  analysed,  using  35  US  Navy  orbits,  and 
fitted  with  theoretical  curves  to  obtain  lumped  values  of  29th-order  harmonic 
coefficients  in  the  geopotential. 
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1 INTRODUCTION 

Ariel  1,  1962-15A,  the  world's  first  international  satellite,  was  launched 
on  26  April  1962  and  decayed  in  the  Earth's  atmosphere  on  24  May  1976  after  a 
lifetime  of  5142  days.  The  satellite  was  cylindrical  in  shape  with  four  paddles 
had  a length  of  0.53m,  a diameter  of  0.58m  and  a mass  of  60kg*.  The  initial 
perigee  and  apogee  heights  were  389km  and  1214km  respectively,  and  the  inclina- 
tion was  53.9°. 

The  orbit  of  Ariel  1 had  an  initial  period  of  100.9  minutes.  Between  1962 

and  1973  the  period  decreased  by  more  than  5 minutes,  under  the  action  of  air 

drag,  and  the  variation  in  inclination  over  these  1 I years  was  analysed  to 

2 

determine  upper-atmosphere  zonal  winds  . During  this  analysis  it  was  noticed 

that  in  November  1969  the  decay  under  the  action  of  air  drag  carried  the  orbit 

through  the  condition  of  29:2  resonance,  when  the  tracks  over  the  Earth  repeat 

every  two  days  after  29  revolutions.  At  the  time  of  this  resonance  a group  of 

2 

35  US  Navy  orbits  revealed  a substantial  perturbation  in  inclination  . 

The  aim  of  this  Report  is  to  evaluate,  for  the  first  time,  harmonic 
coefficients  of  order  29  in  the  geopotential,  from  the  changes  that  occur  in 
the  inclination  and  eccentricity  near  the  time  of  the  29:2  resonance. 

2 THE  THEORETICAL  EQUATIONS  NEAR  29:2  RESONANCE 

The  longitude-dependent  geopotential  at  an  exterior  point  (r,6,A)  may  be 
written  in  normalized  form^  as 

& 

— | P™  (cos  9)  <C„  cos  mA  + S.  sin  mA>N„  , (1) 

r j P 1 Jim  Pm  ) Urn 

P = 2 m= I 

where  r is  the  distance  from  the  Earth's  centre,  0 is  co-latitude,  A is 

longitude  (positive  to  the  east),  p is  the  gravitational  constant  for  the 
3 2 

Earth  (398601km  /s  ),  R is  the  Earth  s equatorial  radius  (6378.1km), 

P?  (cos  0)  is  the  associated  Legendre  function  of  order  m and  degree  P , 

and  C and  S„  are  the  normalized  tesseral  harmonic  coefficients,  of  which 

those  of  order  m = 29  particularly  concern  us  here.  The  normalizing  factor 

N,  is  given  by^ 

Pm  3 

2 _ 2 (2P  + 1)(P  - ml! 

Pm  (P  + m) ! 


(2) 
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The  rate  of  change  of  inclination 

coefficients,  C„  and  S , near  6:n 
’Jim  Jim 


i caused  by  a relevant  pair  of 
. . 4 

resonance  may  be  written 


di 

dt 


— (-t 

sin  1 \a  I 


F„  G„  (k  cos  i - m)<fi 
Jimp  Jlpq 


j£'m+1(Sm"  exPlj^ 


where  F„  is  Allan's  normalized  inclination  function  , G„  is  a function 
Jimp  _ Jlpq  4 

of  eccentricity  e for  which  explicit  forms  are  given  by  Gooding  , denotes 
'real  part  of'  and  j = S-\  . The  resonance  angle  $ is  defined  by  the 
equation 


<!>  = a(u  + M)  + 8 (ft  - v)  , (4) 

where  w is  the  argument  of  perigee,  M the  mean  anomaly,  the  right 

ascension  of  the  node  and  v is  the  sidereal  angle.  The  indices  y,  q,  k and 

p in  equation  (3)  are  integers,  with  y taking  the  values  1,  2,  3 and  q 
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the  values  0,  ±1,  ±2,  ; the  equations  linking  £ , m,  k and  p are  : 

m = y8;  k = ya  - q;  2p  = £ — k . 

Here  6 = 29  and  a = 2 , and  the  m-suffix  of  a relevant  (C^.S^)  pair 

is  given  uniquely  by  the  choice  of  y . The  values  of  Jl  to  be  taken  must  be 
such  that  H » m and  (Jl  - k)  is  even.  The  successive  coefficients  which 
arise  (for  given  y and  q)  may  usefully  be  gathered  together  in  a lumped 
form  and  written  as^ 

q,k  v- » q,k 

C = ) Q C.  , S 

m {_!  xJl  Jim  m 

Jl  Jl 

where  l increases  in  steps  of  2 from  its  minimum  permissible  value  £q  , and 

the  Q.  are  constant  coefficients  with  Q„  = 1 . 

H £ o 

For  the  29:2  resonance  with  y = 1 , we  have  m = 29  and  k = 2 - q ; 

thus  the  affixes  [q,k]  in  equation  (5)  are  [0,2],  [l,l]  and  [- I , 3]  when  q = 0, 

I and  -I  respectively.  Writing  only  the  three  terms  with  (y,q)  * O»0), 

(1,1)  and  (1,-1)  explicitly,  the  theoretical  variation  of  inclination  given 
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by  equation  (3)  may  be  written  for  29:2  resonance  * ’ as 


I 


Q£S£m 


(5) 
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di 

dt 


-JU  (if  [i 

sin  1 \a  / a 


(29  - 2 cos  i)F30>29>14{ 


0,2 


0,2 


S 2 9 sin  $ + C29  cos  3>j 


* 16e(29  - cos 


1,1 


.1,1 


C2g  sin  (4>  - u)  - S2g  cos  ($  - u)| 


+ 126(29  ' 3 C0S  i)f 29,29,131 


-1,3 


.“1,3 


C 2 g sin  ($  + w)  - S2g  cos  (4>  — oj)^ 


q cos  , . . 

+ terms  in  e1’1  . (y4>  - qw) 

sin  n 


) 


(6) 


Only  three  terms  are  given  explicitly  because  it  is  believed  the  others  are 
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small:  terms  with  q * ±2  have  e as  a multiplying  factor  (where  e = 0.04 
for  Ariel  1 at  the  time  of  29:2  resonance),  while  the  terms  with  y = 2 are 
associated  with  harmonics  of  order  58,  which  should  be  much  smaller  than  those 
of  order  29. 


_q,k  _q,k 

The  three  pairs  of  lumped  coefficients  C and  S appearing  in 

m m 

equation  (6)  may  be  written  in  terms  of  the  individual  geopotential  coefficients 
6 

as 


J 0,2 

C29 


= C 


34 


30,29 


32,29,15  M2? 

S \a  / j2,29 

30,29,14  ‘ 30,29,14 


^ (!)V 


29 


(7) 


.1,1 

C29 


C l7F31,29,15  v2 

29,29 


16F 


29,29,14 


l;)V 


18F 


29 


33,29,16 


16F 


29,29,14 


#33. 


29 


, (8) 


and 


.-1,3 

C29 


‘29,29 


13F 


12F 


31,29,14 


29,29,13 


(!)V 


29 


+ 14F33,29, 15  ' " '4 


12F 


29,29,13 


29 


...  , (9) 


and  similarly  for  S , on  replacing  C by  S throughout. 

In  equations  (6)  to  (9)  the  rumerical  values  of  the  three  most  important 
F functions  are: 

on  1/  = 0.247414  sin27i  (15  cos  i-  l)(l  + cos  i)2  = 15.256 x 10  3 for  Ariel  1 
30,29,  14 

— 2ft  — 1 

F,n  „Q  ..  =*  0.506850  sin  1 (1  + cos  i)  = 2.0224  x 10  for  Ariel  1 

F29  29  13  ” 0.443494  sin2^i  (1  + cos  i)3  = 6.8605  x 10  3 for  Ariel  I. 
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be  written 


The  contribution  of  the  (&,m)  harmonic  to  de/dt  for  0:a  resonance  can 
4 


It  = "(a  ) f AmpG£pq  ' ^ ~ i(ke  1 ^ ‘ (Sm  _ 3 V eXH  j (y*  - qu>)}]  . (10) 


For  the  29:2  resonance,  the  theoretical  variation  of  eccentricity  given  by 

_q,k  _q,k  4 5 

equation  (10)  may  be  written  in  terms  of  the  same  (C  , S ) as  ’ 

mm 


s ■ ”(?  j 

■9[  R - 

1 “ a F30,29, 

_0,2 

14e(S29 

sin 

$ + 

-0,2 

C2g  cos  <f) 

16*29,29, 

/c1’1 
14l  29 

sin 

(•  - 

J.«  \ 

(o)  - S29  cos  (4>  - ui)J 

+ I2F29,29, 

ic"1’3 

13129 

sin 

--1,3  > 

+ w)  - S29  cos  (f  + w)> 

+ terms  in  £e 1 q ^ *|q  - J (k  + q)e2  | °°®  (y<f  - qoo)  I . (11; 

Three  terms  are  given  explicitly  in  equation  (11),  those  with  (y,q)  = (1,0), 

(1,1)  and  (1,-1)  . The  main  terms  are  expected  to  be  those  with  (y,q)  = (1,1) 
and  (1,-1)  but  the  (y,q)  = (1,0)  term  is  included  so  that  the  analysis  for 
e conforms  with  that  for  i . 

3 ANALYSIS  OF  THE  29:2  RE SON AM CE 

3. 1 Procedure 

Near  the  time  of  the  29:2  resonance  in  November  1969,  35  sets  of  US  Navy 

elements  were  available  for  analysis.  The  THROE  computer  program  developed  by 
7 4. 

Gooding  ’ provides  a least-squares  fitting  of  equation  (6),  in  integrated  form, 
to  the  observed  variation  in  inclination,  and  similarly  a fitting  of  equation 
(11)  to  the  observed  variation  in  eccentricity,  to  obtain  values  of  the  lumped 

29th-order  harmonic  coefficients  in  the  geopotential.  The  inclination  and 

. 4 

eccentricity  can  be  fitted  simultaneously  using  the  SIMRES  computer  program  . 

3.2  Analysis  of  inclination 

The  35  values  of  inclination,  i , were  cleared  of  zonal-harmonic  and 
lunisolar  perturbations  by  numerical  integration  at  1-day  intervals  using  the 

g 

PROD  computer  program  , and  then  fitted  with  equation  (6)  using  THROE.  A 
standard  deviation  of  0.001°  was  assigned  to  the  35  values  of  i , the  density 
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scale  height  H was  taken  as  58km,  appropriate  to  a height  of  420km  ({H  above 

2 

perigee),  and  the  atmospheric  rotation  rate,  A , was  taken  as  l.lrev/day. 

THROE  can  be  run  for  specified  pairs  of  values  of  (y,q)  and  was  first  run  with 
(y,q)  = (1,0),  (1,1)  and  (1,-1)  , the  terms  given  explicitly  in  equation  (6). 
For  this  run  the  measure  of  fit,  e , had  the  value  1.356.  Several  runs  were 
computed  with  different  (y,q)  terms  added  and  also  with  A changed  to  1.2,  but 
all  these  yielded  higher  values  of  e . So  the  first  choice  of  values  of  (y,q) 
and  A proved  to  be  best.  In  the  initial  run  the  first  three  orbits  and  the 
last  had  values  of  the  resonance  angle  $ more  than  1600°  greater  than  its 
value  at  resonance  (which  was  76°).  So  these  four  outlying  values,  and  one 
badly-fitting  intermediate  value  (at  MJD  40484),  were  omitted:  the  resulting 
fit,  with  30  values,  was  better,  with  e = 1.135  , and  this  run  was  taken  as 
final. 

• # # . , • 

The  values  of  $ and  $ are  shown  in  Fig.l.  The  main  resonance  <1>  = 0 , 

and  the  'subsidiary  resonances'  corresponding  to  the  extra  terms  (y,q)  = (1,1) 

and  (1,-1),  namely  $ = <I>  and  0 = -<L  , all  occur  within  14  days,  between  3 and 

16  November  1969. 

The  30  values  of  i , cleared  of  zonal-harmonic  and  lunisolar  perturba- 
tions, are  plotted  in  Fig. 2 and  the  theoretical  curve  given  by  the  final  THROE 
run  is  shown  with  a broken  line.  The  values  of  the  C and  S coefficients 
obtained  from  the  final  THROE  run  are: 


6-0,2 
10  S29  = 

1.1  ± 0.5 

6-0.2 
10C29  = 

- 0.9 

± 0.4 

6- 1 » 1 

,0C29  " 

- 5.6  ± 3.3 

6-  1 » * 

]0S29  = 

- 6.5 

± 4.5  , 

► 

(12) 

,o6c-‘’3  = 
29 

4.8  ± 1.8 

6--1.3 
10  S29  = 

0.3  ± 

1.5 

3. 3 Analysis  of  eccentricity 

The  35  values  of  eccentricity,  each  with  an  assumed  standard  deviation  of 
0.00004,  were  fitted  with  equation  (11),  in  integrated  form,  using  THROE.  The 
zonal-harmonic  perturbations  are  allowed  for  within  THROE  on  this  occasion;  the 

Q 

lunisolar  perturbations  were  pre-calculated , using  PROD  , but  were  found  to  be 
at  most  0.7  times  the  standard  deviation  and  generally  far  less,  so  they  were 
ignored.  The  density  scale  height,  H , was  taken  as  60km,  appropriate  to  a 
height  of  460km  (3H/2  above  perigee),  and  A was  again  taken  as  l.lrev/day. 


A 
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The  same  pairs  of  values  of  (y,q)  were  used  as  with  the  i fitting,  namely 
(1,0),  (1,1)  and  (1,-1).  The  first  run  revealed  that  there  was  a clear  oscilla- 
tion in  the  US  Navy  values  of  e . So  an  attempt  was  made  to  reduce  the 
oscillation  by  altering  the  value  of  used  in  removing  the  odd  harmonic 

pertubations . By  running  THROE  with  (y,q)  = (0,1)  terms  added,  the  best  value 

of  J was  found  to  be  -3.066  x 10  ^ , instead  of  the  standard  value, 

“6 

-2.531  x 10  . The  use  of  this  value  of  greatly  improved  the  fit, 

reducing  e from  3.61  to  1.78.  A further  run  was  computed,  omitting  the  same 

five  values  as  on  the  final  THROE  run  for  i , for  the  same  reasons.  This 

improved  the  fit,  giving  e = 1.335  . For  all  the  fittings  of  e , the  values 
- - 0 2 

of  (C , S ) 2^  were  expected  to  be  indeterminate  because  they  are  coefficients  of 

small  terms  (see  equation  (11)),  and  only  included  so  that  the  analysis  of  e 
conforms  with  that  of  i . 

The  30  values  of  eccentricity,  cleared  of  zonal-harmonic  perturbations, 
are  plotted  in  Fig. 3 and  the  theoretical  curve  given  by  the  final  THROE  run  for 


e is  drawn  as  a broken 

line . 

The 

values  of  the 

C and 

S coefficients 

obtained  from  the  final 

THROE 

run 

are: 

ft  °’2 

ft-0*2 

10  '29  = 

680  ± 

1080 

10C29  = 

1860  ± 

1010 

ft-1’1 

10C29  = 

- 37  + 

12 

6-  1 ’ 1 

10  S29  = 

- 26  ± 

17 

► (13) 

6-~  1 > 3 
10C29  = 

7.6  ± 

6.6 

6-"1’3 
10  S29  * 

- 13.7  ± 

5.5  . 

J 

The  last  four  values  agree  with  the  last  four  in  equations  (12)  to  within  twice 
the  sum  of  their  standard  deviations,  but  the  first  two  values  in  equations  (13) 
are  obviously  unreliable,  having  standard  deviations  2000  times  greater  than 
the  corresponding  values  in  (12).  An  extra  run  was  computed  without  the  first 
two  terms,  i.e.  with  (y,q)  = (1,0  and  (1,-1),  and  the  values  were  similar  to 
(13),  namely: 

6-  1 > 1 6-  1 ’ 1 
10  C„„  = - 37  + 1 1 10  S„„  = -23+15 


' °6°29  = 

- 37  + 11 

6-  ’ 
10  S29 

6-"1’3 
10C29  = 

9.2  ± 5.7 

6-  ~ 1 ’ 3 
10  S29 

This  confirms  that  the  (y,q)  = (1,0)  terms  had  little  influence  on  the  results. 


r 
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3. 4 Inclination  and  eccentricity  fitted  simultaneously 

The  values  of  i and  e were  fitted  simultaneously  using  the  computer 

4 ... 

program  SIMRES  . With  this  program  there  is  a choice  of  weighting  and  two 

alternatives  were  tried,  first  with  i and  e having  the  same  weight,  and 

second  with  e degraded  by  a factor  equal  to  the  ratio  of  the  final  values  of 

e on  the  THROE  fittings,  namely  1.176  (=  1.335/1.135).  The  second  fitting  gave 

lower  standard  deviations  for  5 of  the  6 constants,  and  is  more  logical  since 

the  pre-assigned  accuracies  in  i and  e were  arbitrary:  the  second  fitting 

was  therefore  preferred.  The  fitting  is  shown  in  Figs. 2 and  3 by  unbroken  lines 

and  the  values  of  the  coefficients  given  by  SIMRES  are: 


A °’2 

A-0’2 

10  S29  = 

1.1  ± 

0.5 

10C29  = 

- 1.0  ± 

0.5 

6- 1 * 1 

I0C29  - 

- 6.6 

± 2.7 

6-’»‘ 

,0S29  = 

- 7.5  ± 

4.6 

► 

(14) 

6--'>3 
10  C29  • 

5.5  ± 

1.7 

6-“ 1 * 3 
>0S29  = 

- 0.2  ± 

1.4  . 

All  these  six  values  agree  with  those  obtained  from  analysis  of  i alone,  as 
given  in  equations  (12),  to  within  about  1/3  of  a standard  deviation,  so  the 
i-fitting  is  dominant.  The  variation  of  e given  by  SIMRES,  as  shown  by  the 
unbroken  line  in  Fig. 3,  differs  considerably  from  the  THROE  fitting,  but  is 
quite  acceptable, 
fitting,  as  Fig. 2 

For  Ariel  1 

_0,2 

C29  = C30,29  " 


_1,1 

C29  = C29,29 


_-l,3 

C29  = C29,29 


The  SIMRES  fitting  of  i is  of  course  similar  to  the  THROE 
shows . 

the  numerical  versions  of  equations  (7),  (8)  and  (9)  are: 


5,0C32,29  + U,5C34,29  ~ 14'0C36,29  + 6,8C38,29 


3,7C40,29  _ 5*4C42,29  ” °'8C44,29  + 3'4C46,29  + 


(15) 


1 1C3 1 , 29  + 44C33,29  ” 92C35,29  + 107C37,29 


47C39,29  “ 36C4 1 , 29  + 4?C43,29  + 9C45,29  + 


(16) 


9C 3 1 ,29  + 26C33,29  ~ 39C35,29  + 23C37,29 


l0C39,29  ~ l9C4 1 , 29  ' 2C43,29  + 13C45,29  + 


(17) 


and  similarly  for  S , on  replacing  C by  S throughout. 
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Equations  (15)  to  (17)  show  that  the  largest  contributions  to  the  lumped 

coefficients  are  likely  to  come  from  individual  coefficients  of  degree  32-43. 

The  expected  order  of  magnitude  of  the  lumped  coefficients  may  be  roughly 

estimated  on  the  assumption  that  the  individual  coefficients  of  degree  Z have 

-5  2 

numerical  values  of  order  10  /l  , so  that  those  with  32  < Z < 43  are  of 
—8 

order  10  . The  numerical  coefficients  in  equation  (16)  are  on  average  about 

eight  times  greater  than  those  in  equation  (15),  while  the  numerical  coeffici- 
ents in  equation  (17)  are  about  three  times  greater  than  those  in  (15),  so  the 
magnitudes  of  the  three  successive  lumped  coefficients  would  be  expected  to  be 

in  the  ratio  1:8:3.  The  actual  ratios  from  equations  (14)  are  <J(  1.1^  + 1.0^): 

V2  2 / 2 2 

(6.6  + 7.5  ):  v(5.5  + 0.2~)  or  1:7:4  approximately;  so  the  agreement  is 

excellent. 

From  equations  (15)  to  (17)  the  expected  orders  of  magnitude  of  the 
numerical  values  of  the  three  lumped  coefficients  are  0.2  x 10  1.5  * 10  ^ 

and  0.5  * 10  ^ , The  actual  values  obtained  in  equations  (14)  are  within  a 
factor  of  5 of  the  values  indicated  by  this  rough  order-of-magnitude  estimate, 
the  actual  values  being  generally  the  greater. 

4 CONCLUSIONS 

The  perturbation  noticed  in  the  orbital  inclination  of  Ariel  1 near  29:2 
resonance  has  been  analysed,  together  with  the  perturbation  in  eccentricity,  to 
obtain,  for  the  first  time,  values  of  lumped  29th-order  harmonic  coefficients 
in  the  geopotential.  The  values  are  given  in  equations  (14),  the  symbols  being 
defined  in  equations  (1)  to  (9).  The  numerical  values  derived  are  not  particu- 
larly accurate,  but  are  as  good  as  can  be  expected  in  view  of  the  limited 
accuracy  of  the  orbital  data. 

The  success  of  this  analysis  shows  that  accurate  values  of  29th-order 
harmonics  should  be  obtained  in  future,  if  satellites  which  pass  through  29:2 
resonance  are  intensively  observed  so  that  accurate  orbits  can  be  determined. 
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